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Abstract 

We apply the inverse scattering method to the midi-superspace models that are char- 
acterized by a two-parameter Abelian group of motions with two spacelike Killing vec- 
tors. We present a formulation that simplifies the construction of the soliton solutions 
of Belinskii and Zakharov. Furthermore, it enables us to obtain the zero curvature 
formulation for these models. Using this, and imposing periodic boundary conditions 
corresponding to the Gowdy models when the spatial topology is a three torus T'^, 
we show that the equation of motion for the monodromy matrix is an evolution equa- 
tion of the Heisenberg type. Consequently, the eigenvalues of the monodromy matrix 
are the generating functionals for the integrals of motion. Furthermore, we utilise a 
suitable formulation of the transition matrix to obtain explicit expressions for the in- 
tegrals of motion. This involves recursion relations which arise in solving an equation 
of Riccati type. In the case when the two Killing vectors are hypersurface orthogonal 
the integrals of motion have a particularly simple form. 



1. Introduction 



The Einstein field equations for space-times admitting a two-dimensional Abelian 
group of isometries acting orthogonally and transitively on non-null orbits are non- 
linear partial differential equations in two variables. For timelike orbits the equations 
are elliptic, whereas for spacelike orbits the equations are hyperbolic [1] . Although the 
space-times admitting two commuting Killing vectors are not the most general, they 
can represent interesting physical situations with stationary axial symmetry, planar 
symmetry or cylindrical symmetry [2]. 

Since the pioneering work by Geroch [3], it has been known that the field equations 
in the stationary axisymmetric case admit an infinite dimensional group of symmetry 
transformations. This result has encouraged research in solution-generating methods, 
the idea being that the complete class of solutions can be generated from a particular 
solution, such as flat space [4]. Subsequently, several solution-generating techniques 
have been developed, such as the Kinnersley-Chitre transformations [5], the Hauser- 
Ernst formalism [6], Harrison's Backlund transformations [7] and the Belinskii and 
Zakharov inverse scattering method [8]. The relations between the different approaches 
were discussed by Cosgrove [4], with Kitchingham [9] adapting these methods to the 
hyperbolic case. 

Our paper uses the framework of the Belinskii and Zakharov inverse scattering 
method, for the case when the field equations are hyperbolic. The inverse scattering 
method is a powerful method for solving certain systems of non-linear partial differ- 
ential equations. The main step in this formalism is to write down linear eigenvalue 
equations whose integrability conditions are the given non-linear system. The methods 
of functional analysis can be applied to generate new solutions of the linear system 
from old, and hence new solutions of the original system from old [10]. The particular 
solutions that can be generated are the soliton solutions. The soliton solutions share 
a number of common properties with classical particles, namely, they are localized 
solutions that propagate energy, have a particular velocity of propagation and some 
persistence of shape which is maintained even when two solitons collide [2]. As shown 
by Belinskii and Zakharov, the soliton transformation needs to be generalized when 
applied in the context of the Einstein equations with a two-parameter Abelian group 
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of motions. The generalization is that the stationary poles are substituted by the 
pole trajectories. We will present an equivalent formulation of the linearized system. 
Similarly to Belinskii and Zakharov [8], our linearized system is defined with the use 
of two differential operators which involve derivatives with respect to the complex pa- 
rameter A. We then define a new complex parameter a;, and show that this simplifies 
the linearized system. We also point out the properties of the map between the two 
complex parameters u>{t, z, A) and X{t, z, ou). We then show that our formulation yields 
soliton solutions equivalent to the original solitons of Belinskii and Zakharov. 

We next construct the zero curvature formulation for the system. The zero cur- 
vature formulation is an important characteristic of integrable systems. A direct con- 
sequence of the zero curvature formulation, for a given system, is the fact that, when 
periodic boundary conditions are imposed, the equation of motion for the monodromy 
matrix is an evolution equation of Heisenberg type [11]. Hence, the eigenvalues of the 
monodromy matrix are conserved. In the context of the two spacelike Killing vector 
reduction of general relativity, periodic boundary conditions amount to the compact- 
ification of the z direction into a circle. This then corresponds to the Gowdy models 
when the spatial topology is a three torus [12]. 

The next step is to obtain a suitable parameterization for the transition matrix. 
To achieve this we have to solve a system of four partial differential equations. This 
problem reduces to one of solving two equations of Riccati type. We then look for 
solutions to the Riccati equations in the form of power series in (Ail). The solutions 
are given through recursion relations. For integrable systems with fixed poles, the 
integrals of motion are then given by the coefficients in the Laurent expansion of the 
eigenvalues of the monodromy matrix around the poles (Ail). However, in the case 
of the Einstein field equations for space-times admitting a two-dimensional Abelian 
group of isometries acting orthogonally and transitively, as we have remarked earlier, 
the fixed poles are substituted by the pole trajectories, i.e., {X{uj,t,z) i 1 ). Due to 
the fact that A is time-dependent, i.e., dtXico, t, z) ^ 0, we are not able to identify the 
"local" integrals of motion as the coefficients in the expansion of the eigenvalues of 
the transition matrix. Instead, "local" integrals of motion are given as the eigenvalues 
of the transition matrix for fixed values of the complex parameter u in the domain in 
which all the relevant algebraic series converge uniformly. In the case when the Killing 
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vectors are hypersurface orthogonal the integrals of motion have a particularly simple 
form. 

This paper is organized as follows. In section two we formulate the inverse scat- 
tering method as applied to the two Killing vector reduction of general relativity. This 
involves the feature that the derivatives defining the first-order formulation of the 
equations of motion also involve derivatives with respect to the spectral parameter. 
We show how this may be dealt with by defining a new complex parameter, and we 
discuss some properties of this map. We then show how the approach of Belinskii and 
Zakharov may be adapted to our case, and thus we obtain the soliton solutions. In 
section three, we turn to the zero curvature representation of the equations of motion, 
using this to define a transition matrix in the usual way. The integrals of motion are 
then seen to be given in terms of the related monodromy matrix. In section four, we 
utilise a suitable formulation of the transition matrix to obtain explicit expressions 
for the integrals of motion. This involves recursion relations which arise in solving an 
equation of Riccati type. Finally, in section five we present our conclusions. 

2. The Inverse Scattering Method and the Soliton Solutions 

We will consider the midi-superspace models that are characterized by the existence 
of a two-parameter Abelian group of motions with two spacelike Killing vectors (the 
case when one Killing vector is timelike and the other spacelike is similar and we will 
not consider it separately). Let us choose coordinates adapted to the action of the 
symmetry group so that the metric assumes the following form [8] 

ds'^ = -fdt^ + fdz^ + gabdx^dx^ (2.1) 

where a,b = 1,2, {x'^ , , , x^} = {t,x,y,z}, / is a positive function and gab is a 
symmetric two-by-two matrix. The function / and the matrix gab depend only on the 
co-ordinates {t, z}, or equivalently on the null co-ordinates {^, ry} = {^{z+t), ^{z — t)}. 
There is a freedom to perform the co-ordinate transformations 

{c,v}^m),m}- (2.2) 

It is easy to see that the transformations (2.2) preserve both the conformally flat 
two-metric f{—dt^ + dz"^) and the positivity of the function / if d^^ drjfj > 0. 
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The complete set of vacuum Einstein equations for the metric (2.1) decomposes 
into two groups of equations [8]. The first group determines the matrix gab and can 
be written as a single matrix equation 



d^{ad^gg-') + d^{ad^g g-') = , (2.3) 

where = det g and 77} are the null co-ordinates. The second group of equations 
determines the function /(^, rj) in terms of a given solution of (2.3): 

9? (In a) 1 „ 
d^^^f) = /tt^ + 7 TrA' , (2.4a) 

d?,(lna) 1 „ 
d,{lnf) = ^j-^ + TrB^ , (2.46) 

where — d^a, — drja and the matrices A and B are defined by 

A = -ad^g B = adr^g g~'^ . (2.5) 

The dynamics of the system is thus essentially determined by eqn. (2.3) and for this 
reason we will concentrate on it in the following. 

Taking the trace of eqn. (2.3) and using the definition for a, we obtain 

a^rj = Q- (2.6) 

The two independent solutions of this equation are 

a = aiO + h{Ti) , (2.7a) 
(3 = a(0 - 6(r/) . (2.76) 

Using the transformations (2.2), one can bring the functions a(^) and 6(77) to a pre- 
scribed form. However, we will consider the general form without specifying the func- 
tions a(^) and 6(77) in advance. 

Let us now consider the system of equations 

Vr,A + V^B = 0, (2.8a) 
dr,A - d^B + [A, 5] = , (2.86) 
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where = + (x^a~^, = 5^ + a^o;"^ and [ , ] denotes the commutator in 
the Lie algebra of the group GL(2,R). The general solution of the equation (2.8b) is 
of the form 

A = 8^11-'^, B = dr,ir\ (2.9) 

where I is an element of the group GL{2,R). In addition we impose the constraint 
= detL Eqn. (2.6) follows from the trace of eqn. (2.8a), once we substitute eqn. 
(2.9) into (2.8a). However, we still have more degrees of freedom in eqn. (2.8a) than 
in (2.3). We thus need to impose some additional conditions in order to recover the 
correct number of degrees of freedom. For this reason we impose the constraint I = T , 
where l'^ is the transpose of the matrix /. It then follows, once the additional constraints 
are imposed, that eqns. (2.8) are equivalent to eqn. (2.3). Due to the constraint I = l'^, 
the matrices A, B in eqn. (2.9) must be taken to be only those which can be written 
in this form using a symmetric matrix I. The sets of such matrices A, B form subsets 
(but not subgroups) of GL{2,R). 

The crucial step in the inverse scattering method is to define the linearized system 
whose integrability conditions are the equations of interest, in our case eqns. (2.8). 
Following ref. [8] , we first define the two differential operators 

D^ = d^ -^^Aa;, , (2.10a) 
D.^d^ -^^A^A , (2.106) 

where A is a complex parameter independent of the co-ordinates {C,r]}. Notice that 

and Drj are invariant under the coordinate transformation {A, ^,7/} — > {jj^jV}- 
is straightforward to see that the differential operators and Dr/ commute since a 
satisfies the wave equation (2.6) 

P..^5,l = ^(^-^)A9. = 0. (2.11) 
The next step is to consider the following linear system 

D^^=-J^i^^ (2.12a) 

^'?^=aTi'^' ^^'^^^^ 
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where il^{X, ^, rj) is a complex matrix function, and the real matrices A, B and the real 
function a do not depend on the complex parameter A. The integrability conditions 
for the system (2.12) are eqns. (2.8). To prove this we first apply the operator 
to eqn. (2.12a) and subtract from this the result of applying the operator to eqn. 
(2.12b). The left-hand side vanishes using eqn. (2.11) and the right-hand side is a 
rational function of A which vanishes if eqns. (2.8) are satisfied. In order to take into 
account the additional constraints that we have imposed we require that 

V;(A)='0(A), (2.13a) 
g = ^{^)r{\), (2.136) 

where A is the complex conjugate to A, g{^,f]) is a symmetric two-by-two matrix of 
functions, whose determinant is a^, and ip{X) satisfies eqns. (2.12). The condition 
(2.13b) implies that A and B have the correct number of degrees of freedom. To see 
this, we apply the operator to equation (2.13b). Using the fact that ip satisfies eqn. 
(2.12a) we obtain 

d^g=j^(^XAg-gA^y (2.14a) 

Similarly, using (2.12b), we obtain 

dr,g-^j^{>^Bg + gB-y (2.146) 

Taking the transposes of the right-hand sides of equations (2.14), we deduce that 
Ag = gA'^ and Bg = gB'^. Consequently, A and B have the form 



A = d^gg-\ B = d^gg-\ (2.15) 

The standard application of the inverse scattering method to field theories in 
(l-l-l)-dimensions differs from the present situation in that the linearized system does 
not usually involve differentiation with respect to the spectral parameter A, as occurs 
here. The present situation can be improved by defining a new complex parameter u) 

by 

u=\[j + 2(5 + aX). (2.16) 
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A straightforward calculation shows that D^u = D^oj = 0. Furthermore, by making 
a co-ordinate transformation {A, ^, 77} {u, ^, 77} we reduce and to and dr^ 
respectively. In order to write the linearized system (2.12) in the new co-ordinates we 
have to invert the relation (2.16), i.e., we have to know A as a function of {uj,^,t]). 
The inverse of (2.16) is not unique - we encounter the following two solutions for A: 

^-(^)-((^)^-0*- 

At this point, we have a choice of two different co-ordinate systems, the first corre- 
sponding to A_|_ and the second to A_. We can now write the linearized system in the 
new co-ordinates as 

9i^ = -J^i^^ (2.18a) 
^'?^ = T^^' (2.186) 

where 'ip{(jj,^,r]) and A± is given by (2.17). It is straightforward to check that the 
integrability conditions for this system are equations (2.8). To show this it is necessary 
to use equations (2.7) and (2.17). 

Let us point out some properties of the map (2.16). The transformation A ^ 
leaves (2.16) invariant. The A plane is mapped into the two-sheeted Riemann surface 
which covers the entire A plane with the branch points at u = P ± a. The map (2.16) 
takes the circle |A| = p, into the ellipse Cp given by the parametric equation 

a 1 

w= -(p+ -)cos(^ + /3, (2.19a) 
2 p 

a 1 

v=-{p--)smcf>, (2.196) 

where u = u + iv and is the phase of A, i.e., A = pe^'^. In particular, the image of 
the circle |A| = 1 is a closed interval on the real axis in the u; plane. From our previous 
discussion it follows that the two formulations, the first defined by (A, ^,77) and the 
linearized system (2.12), and the second defined by (a;,^,?7) and the linearized system 
(2.18) are completely equivalent. In what follows we will use both formulations. 

It is straightforward to obtain the soliton solutions following the inverse scattering 
method [8]. We begin by setting ip\^^^^ = I. Then, from eqn. (2.13) it follows that 
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S'(C)^) = i^i^iCiV)- III the inverse scattering approach we start with a given solution 
fl'o(C)^) of eqns. (2.3). From eqn. (2.15) we then determine ^0(^,77) and Bo{^^r]). 
Substituting Ao{^,ri) and Bo{$^,r]) into eqn. (2.12) and solving the linearized system 
we obtain ipo{X, ^, rj), with go{^, t]) — ^o(0, C,t]) ^ ■ We now define x cts 

^{\^,v) = xi\^,v)M\tv)- (2.20) 
The linearized eqns. for from eqns. (2.12), are 

We also have the constraint on x, from eqns. (2.13), 

1 

f (C, ^) = X(^, V) 9o{^, V) X^{\ ^, V) 
Setting x\x=±oo = ^ (2-22) it follows that 

^(C, ^) = X(0, e, ^) yo(C, ^) • (2.23) 

The soliton solutions arc characterized by the points in the A plane at which the 
determinant of x is equal to zero and x~^ has a simple pole, and similarly, the points 
at which the determinant of x~^ is equal to zero and x has a simple pole [8]. Thus, x 
and x~^ cire rational matrix functions of A with a finite number of simple poles 

X 

where the matrices Tit and Sk do not depend on A, and fik is the complex conjugate of 
Hk, and Uk of Uk- From eqn. (2.24a), it is easy to check that the condition x\x-±oo ~ ^ 
is satisfied. Also, from eqn. (2.22) it follows that the positions of the poles of x and 

^ Integration of the linearized system is straightforward in the case when the metric 
fl'o(C) v) is diagonal [13]. However, when go{$,, rf) is non-diagonal this step is non-trivial, 
an example of this being the case of the Bianchi type II models [14,15]. 



(2.21a) 
(2.216) 

(2.22) 



^ c c 
^-^ W — Uh X — Vhf 



(2.24a) 

(2.246) 
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the poles of x -"^ are related by //fc^^fc = 1) k = 1, ....,N. This implies, using eqn. (2.16), 
thata;fe = i(^ + 2/3 + a//fc) = i + 2/3 + az/fe) , k = l,....,N. 

The next step is to rewrite eqns. (2.21) in the more convenient form 

= {Da)x-' - X^X-' , (2.25a) 



A + 



Y = (^.x)x"'+X^X-'. (2.256) 



Substituting eqn. (2.24) into (2.25) and setting the residuals at the poles A = /Xfc of 
the right-hand-side equal to zero, we obtain a set of N first order equations for the 
matrices TZ^, k = 1, ....,N. Similarly, from the identity xx""*^ = I we obtain the system 
of N algebraic equations 

TZkX'HlJ'k) = , (2.26a) 
and from eqn. (2.22) we obtain another system of N algebraic equations 

7^fe^o(I + E(^^ + ^^)) =0, (2.266) 

where we have used the identity jJ.ki'k = 1. We also have similar equations for the 
matrices TZk and the poles A = Jlk-, k = 1,...,7V. This yields the n-soliton solutions 
equivalent to those of Zakharov and Belinskii, [8]. In addition, equations (2.4), which 
determine the function /, can then be integrated explicitly following ref. [8,16]. We 
will not present these results here. Instead, we will turn to the study of the zero 
curvature formulation and the integrals of motion. 

3. The Zero Curvature Formulation and The Integrals of Motion 

We now construct the zero curvature formulation for equations (2.8) and, using the 
techniques of ref. [11] in the case when the z direction is compactified into a circle S^, 
obtain the generating functional for the integrals of motion. Our first step is to perform 
a co-ordinate transformation {^, ry} — > {t, z} and define U{t, z, A) and V{t, z, A) by 

^('■^■^) = 5(-A4T + AfT)- 

''('•^•^)=i(-A^-AfT)- 
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The equation of motion for U and V is 



DtU - D,V +[U,V]=0, 



(3.2) 



where A = | - Drjj , = \ [D^ + Dr,j , with and given by eqn. (2.10). 
Substituting eqn. (3.1) into (3.2), it is straightforward to show that eqn. (3.2) is 
equivalent to eqns. (2.8). 

The fields U{t,z,u) and V{t,z,uj) are expressed, in the co-ordinates {t,z,uj}, by 
the formulae that can be obtained from (3.1) by the substitution A — > X+{t, z,uj), 
where X+{t,z,u>) is defined in (2.17) ^ . Consequently, the equation of motion in this 
co-ordinate system is the zero curvature equation, i.e., eqn. (3.2) with the substitution 
Df — > dt and Dg ^ d^. We will consider the case when periodic boundary conditions 
are imposed 



Notice that we impose periodic boundary conditions on all fields, including a and f3. 
These boundary conditions correspond to compactifying the z direction into a circle 
S^. For example, this corresponds to the Gowdy models when the spatial topology is 
a three torus [12]. We consider the transition matrix 



U{t,z + 2L,u) = U{t,z,uj) , 
V{t,z + 2L,oj) = V{t,z,Lj). 



(3.3a) 



(3.36) 




(3.4) 



The transition matrix satisfies 



dz^T{t,z-i_,zo,u)) = U{t,zi,uj)T{t,zi,zo,u) 



(3.5) 



with the condition 



T{t,zi,zo,uj) 



= I. 



(3.6) 



Zo=Zl 



We apply dt to eqn. (3.5), obtaining 



dtid,T)=dtiU)T+UdtiT). 
^ From hereon we will restrict ourselves to the plus sign in eqn. (2.17). 
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(3.7) 



Using the equations of motion for {U, V}, the zero curvature equation (3.7) can be 
written 

{dtT - VT) = U{dtT - VT) . (3.8) 

It foUows from (3.5) that 

dtT{t,zi,zo,uj) = V{t, zi,u>)T{t,zi,zo,Lu) 

+ T{t,zi,zo,uj)C{t,zo,oj), (3.9) 

and, iising the condition (3.6), we get C{t, zqjUj) = —V{t, zo,u>). Thus the equation of 
motion for the transition matrix is 

dtT{t,zi,zo,u)) = V{t,zi,uj)T{t,zi,zo,u}) 

-T{t,zuZo,uj)V{t,zo,uj), (3.10) 



We now define the monodromy matrix to be the transition matrix afong the 
fundamental domain —L < z < L, i.e., 

TL{t,oj) = T{t,L,-L,u;). (3.11) 

From eqn. (3.10), using the periodic boundary conditions (3.3), it foUows that the 
equation of motion for the monodromy matrix TL{t,(ju) is an evolution equation of 
Heisenberg type 

dtTL{t,uj) = [V{t,L,u;),TL{t,<.v)] . (3.12) 

This implies that the eigenvalues of the monodromy matrix T£,(i, cj) are conserved, or 
equivalently 

dtTrTL{t,uj) = 0, (3.13a) 
dtTT{TL)\t,uj) = 0. (3.136) 

Our conclusion is that the functions 

EL{uj) = TTTL{t,u), (3.14a) 
FL{u) = TT{TL)\t,uj), (3.146) 
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are the generating functions for the conservation laws. This result is a direct con- 
sequence of the zero curvature formulation and the periodic boundary conditions as 
shown above. Similar results are known for most if not all of the dynamical systems 
that admit a zero curvature formulation [11]. The standard method to obtain the ex- 
plicit expressions for the integrals of motion involves solving the equations of Riccati 
type. The integrals of motion are then identified as the coefficients in the Laurent ex- 
pansions of the generating functions. However, as we will show in the next section, in 
our case the generating functions (3.14) depend on the complex parameter ui through 
A-).(t, z, (jj). Consequently, we cannot identify the integrals of motion as the coefficients 
in the expansion of the generating functions. Instead, our integrals of motion are given 
by algebraic series, for those values of uj for which the series converge uniformly. 

4. The Integrals of Motion 

We can now obtain explicit expressions for the integrals of motion. We first write the 
transition matrix in the form 

T(t, ^1 , ^0, a;) = (l + W{t, z^.u)) e^^*'^^'^"''^) (l + W{t, zq, o;)) , (4.1) 

where I is the identity two-by-two matrix, W{t, z, u) is an off-diagonal matrix and 
Z{t, zi, zq.lj) is a diagonal matrix [11]. Substituting eqn. (4.1) into (3.5) we obtain 
the following system of equations 

d^^W^^) + W^(3) (^a.,Z(i) - c»,,Z(2)) - + C/(2))t^(3) _ [/(3) = , (4.2c) 

+ ^(4) (^Q^^z^l) + Q^^z^2)^ - - C/(2)^^(4) _ ^(4) ^ Q _ ^4_2cZ) 

Here we have used the notation Z = Z(^) r(i) + Z^^) r(2), W = W^^^t^^) + 1^(^V(4) and 
U = Yl'i=i U^^^ ^i)- choice for a basis in the Lie algebra of the group GL(2,R) is 

T(l) = I, T(2) = il -l)' ^(3) = (S i) and T(4) = (; S)- 

Substituting eqns. (4.2a,b) into eqns. (4.2c,d), we obtain the system of equations 
that determines W{t,zi,uj) 

^ ^(4) (w^3)y _ 2[/(2)^(3) _ [/(3) = , (4.3a) 
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+ + 2^7(2)^(4) _c/(4) = o, (4.36) 

with periodic boundary conditions on W{t, zi,uj), i.e., W{t, zi,uj) = W{t, z\ + 2L, a;). 
Once we have solved eqns. (4.3) for Wit^zx^uj)., eqns. (4.2a,b), together with the 
boundary condition Z{t, Zi, zo,uj)\^^_^^ = 0, determine Z{t, zi, zq^co): 

Z^'\t,z,,zo,uj) = Hdzlu^'^ + Vc/(3)W^W + ^ ^ (4 4^) 

Z^^\t,z,,zo,u;) = rdz(u^^^ + l(^c/(3)p^(4) _ [/(4)|^(3)^ ^ _ (4_4^) 

The main problem is thus to obtain the solutions to the system (4.3). These 
equations are of Riccati type. Given an equation of Riccati type with arbitrary coef- 
ficients together with a particular solution, it is possible to reduce the equation to a 
linear first order system [17]. However, we do not have particular solutions to equa- 
tions (4.3), so we need a different approach. We make a coordinate transformation 
{t, zijUj} — > {t, zi, A} and as a result we obtain the equations (4.3) in the form 

D,,W^^^ + C/W (W^^^y - 2U^^^W^'^ - C/(^) = , (4.5a) 
i?,, + [/(3) (w^^)y + 2[/(2)t^(4) _ c/(4) ^ , (4.56) 

We now expand the fields W^^'^ and W^'^'^ as power series in (A — 1) 

oo oo 

M/(3) = Wi"^^ (A - 1)'^ , W^^^ = Wi""^ (A - 1)'^ ■ (4.6) 

n=0 n=0 

Substituting the first equation of (4.6) and (3.1a) into (4.5a) and using the expansion 
= |E^q(— |)"(A — 1)" we obtain the recursion relation (TV = 0, 1, 2, . . .) 

(N + 2) ^wSU = a.. - ^ E wi^^w'^l_„ + ^ f; (- i) " X 

n=0 m=0 

N—m „ AT" 1 1 

E i^-f '<^l.-» - (f + - ^<^')<+i - f f - J E (4) X 

n=0 n=0 

n=l 
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together with 

^W<« = -^(H/<'')V^«VFr' + ^. (4.8) 

Notice that Wq^^ is arbitrary, and for every choice of Wq^^ we have a different solution. 
From eqn. (4.5b), we obtain similar equations determining the Wn^\ which may be 
obtained from eqns. (4.7) and (4.8) by the replacements W"i^^ ^ ^ A^^\ 

5(4) ^ 5(3)^ ^(2) ^ _^(2) 5(2) ^ _5(2)_ 

In this way we have obtained the solutions to equations (4.5) in an open neigh- 
borhood of the point A = 1. Our next step is to perform a coordinate transformation 
{t, zi,X} ^ {t, zi,(jj}. Then W = W{t, z, X+{t, z, oj)) becomes a function of t, z and cu, 
and eqn. (4.6) becomes an algebraic series in powers of (A+(t, z,u) — 1). Substituting 
this expansion and the expansion j^j^ = \ Xl^o "1)"^ i^^to (4.4) we obtain 

J Zq + 

n=0 

+ 2 E {-IT + B^H-'IU) )(A, - D") .(4.9«) 

m=0 

+ i E (^''' (4) - (■*'"<'' - ■*""'»'"+0 

n=0 

+ 5 E (-i) (B'^W^iU - SW^-iiU) )(A+ - D") .(4.96) 

m=0 

where wi^^ and M^i^^ are determined from eqns. (4.7) and (4.8). Thus we have 
obtained the functions Z^^\t, zi, zo^u!): defined on open neighborhoods of the hyper- 
surface X+{t,z,u!) = 1 where the relevant algebraic series converge uniformly. 

We now proceed to construct the functions Z^*) {t, zi,zo,uj), defined on open neigh- 
borhoods of the hypersurface X+{t, z,uj) = —1. Our first step is to construct the solu- 
tion to the equation (4.5) in an open neighborhood of the point A = — 1. We expand 
the fields W^^'^ and W^^'^ in powers of (A -|- 1) 

oo oo 

= ^ Wf)(A+l)^, W^(^) = J]VF^^)(A+1)", (4.10) 

n=0 n=0 
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and the expansion = — |X]^o(^2^)" l^ads to the following recursion relation 
(JV = 0,1,2,...) 



n=0 m=0 
N — m N 

n=0 n=0 



together with 

^t^f ) = ^(W^'^)' - B^^)W^'^ - \b^') . (4.12) 

As with eqn. (4.8), is arbitrary, and for different choices of we have different 
solutions. From (4.5b) and (4.10) we can obtain similar equations determining the 
wlt\ by the following replacements in those equations: Wn^ — > wlf^, A^^^ <-> A^^\ 

5(4) ^ 5(3)^ ^(2) ^ _^(2) 5(2) ^ _5(2)_ 

We now perform the coordinate transformation {t, A} — > {t, ^, a;}. Then (4.10) 
becomes an algebraic series in powers of (A+ + 1). Using the expansions (4.10,11) and 
= -\ T.n=Q{^T we obtain 



zW(^,.„.o,u;) = jy^C^X^A''^'^ + \{b^M'^ + B^^W^'^) ) 
+ I E (-^'''f + [B'^'wltl + B''>WS>,) 

n=0 

+ E (^^'^^n-- + Ai^'wi'}^) ) (A+ + 1)-) , (4.13a) 

m=0 
n=0 

+ E ^ (^^'^^iiU - A^'^wi'lJi ) (A+ + 1)-) , (4.136) 



m=0 
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where Wn^^ and Wn^^ are determined from the recursion relations derived for them 
above. We have thus determined the functions Z^'^\t, ^i, 2:0,0;) in an open neighbour- 
hood of the hypersurface A+(t, z, uj) = —1 where the relevant algebraic series converge 
uniformly. 

Since we know (locally) the functions Z^'^\t, zi, zq^uj), we can obtain expressions 
for the integrals of motion of our system. Using (4.1) and the fact that W{t,z^(jj) = 
W{t, z + 2L, (jj) it follows that 

El{uj) = TrTL{t,uj) = Tre^^-^ '-•^''^^ , (4.14a) 

From (4.14), after a straightforward calculation, it follows that 

zW(L,-L,.)=il„(^^!M_ZiM), (4.15a) 

ZP)(L,-L,..)=icosh-'(-^/iM_^). (4.156) 

2 \El (cu) - Fl{oj)^ 

Consequently, from (3.13) and (4.15), we have 

dtZ^'\L,-L,uj) = 0. (4.16) 

The conclusion is that the Z^^\L, — L,a;) are the integrals of motion for every fixed 
value of UJ which belongs to the domain in which the relevant series converge uniformly. 

Our final remark is on the case when the two Killing vectors are hypersurface 
orthogonal. In our formulation, this case corresponds to the vanishing of the off- 
diagonal matrix W in (4.1). The transition matrix T is then diagonal. Consequently, 
our integrals of motion have a particularly simple form 

Z^^\L,-L,lj) = J dzU^^\ (4.17a) 



with 



Z^^\L,-L,u) = J dzU^^\ (4.176) 
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To obtain (4.18) we have used equations (2.15) together with the parametrization 

g = e«^(i)+Tn2) . (4.19) 

Work on the geometrical interpretation of the integrals of motion Z'^^\L^ —L,u) 
and the boundary conditions that correspond to gravitational plane waves and gravi- 
tational cylindrical waves is in progress [18]. 

5. Conclusions 

In this paper we have formulated the inverse scattering method as applied to the midi- 
superspace models characterized by a two-parameter Abelian group of motions with 
two spacelike Killing vectors. The application of the inverse scattering method to this 
model involved the feature that the derivatives defining the first-order formulation of 
the equations of motion also involved derivatives with respect to the spectral parameter 
A. We dealt with this by introducing a new spectral parameter uj. We also discussed 
the properties of the map between the two complex parameters u}{t, z, A) and A(t, u)). 
Then we showed how the approach of Belinskii and Zakharov could be adapted to our 
case, and hence found the soliton solutions. 

Our next step was to obtain the zero curvature representation of the Einstein field 
equations for space-times admitting a two-dimensional Abelian group of isometries, in 
the case when the orbits are spacelike. Then, using periodic boundary conditions which 
in the present context correspond to the three torus Gowdy models, we showed how 
the zero curvature formulation implies that the equations of motion for the transition 
matrix are of Heisenberg type. Consequently the eigenvalues of the transition matrix 
are conserved. 

To obtain explicit expressions for the integrals of motion we had to solve a system 
of four partial differential equations. This problem was reduced to the problem of 
solving two equations of Riccati type. The solutions of these equations were given 
through recursion relations. The final results are the integrals of motion which are 
given as powers series in ( A_|_(t, z,^;) ± 1 ). In the case when the two Killing vectors 
are hypersurface orthogonal the integrals of motion have a particularly simple form. 
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Future research will include the elaboration of the geometrical interpretation of 
the integrals of motion that we have obtained, as well as an investigation of different 
boundary conditions, such as those corresponding to planar symmetry and cylindrical 
symmetry [18]. The problem of quantization of the soliton solutions on the Bianchi 
type II background is under investigation and we will report our results shortly [15]. 
Finally, our main goal is to obtain the quantum theory for the two Killing vector 
reduction of general relativity [19]. 
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